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Introduction 



j> ■ Since the 1970s, the mathematics of finance has been developing in the direction of creating progres- 
OO sively more general models for price evolutions of basic financial assets, from Brownian motion to the 
' almost all-encompassing semimartingale processes. This has led to the creation of theory of no-arbitrage 
^ . and completeness, as well as hedging and optimal investment (see Merton Duffie Delbaen and 
Schachermayer jH, Karatzas and Shreve [l^, etc). Although it would be difficult to improve these theo- 
' retical findings in terms of structural generality, the efforts to calculate exact theoretically and practically 
C , significant formulas for option pricing and corresponding optimal investment have been successful only for 
' those models of financial markets in which the increments of underlying random processes are independent 
I> (Wiener, Poisson, Levy processes etc.). 

^ Meanwhile, the development of non-semimartingale models focused mainly on accounting for the depen- 
■ dence of asset prices on the past (long-term memory processes, fractional Brownian motion etc.). However, 
as of today there is still no commonly accepted theory on this topic, nor adequate uses of existing theoretical 
results in practice (see, for instance, Bjork and Hult 0]). 

Another direction, which lies somewhat in between the two trends mentioned above, utilizes markovian 
dependence on the past and the technique of Markov random processes. In our opinion, this approach 
and the corresponding methods (see Elliott and van der Hoek [l3]) are still inadequately reflected in 
contemporary flnancial mathematics. 

Our paper deals mainly with this direction of study of financial markets. As a basis for building the 
model, we take a Markov process a{t) with values ±1 and transition probability intensities \±. Using these, 
we define processes Co-(() = c-t, h^ii-^ = h± and r^(j) = t±, r-t > 0. Let us introduce X'^{t) = Jq c„(^g^ds and 
a pure jump process J" = J"{t) with alternating jumps of sizes h±. The evolution of the risky asset S{t) 
is determined by the stochastic exponent of the sum -|- J*^. The risk-free asset is given by the usual 
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exponent of the process = ^^{t) = r^^^g-^ds. Here and below the superscript a indicates the starting 
value a = (t(0) of ait). 

In view of such trajectories, the market is set up as a continuous process that evolves with velocity c+ 
or c_, changes the direction of movement from c± to c=p and exhibits jumps of size h± whenever velocity 
changes. 

The interest rate in the market is stochastic with the values r± such that {c± — r±)h± < which means 
that the current trend of discounted prices and the direction of the next price jump should be opposite. 
A process X'^ = X'^{t) is defined by a pair of states (c±, X±) and is called a telegraph process with states 
(c±) It describes continuous price trends (upward or downward) between random instants. Changes 
in these trends are accompanied by jumps of sizes h±. 

Our model uses parameters c± to capture bullish and bearish trends in a market evolution, and values 
h± to describe sizes of possible crashes, jumps and spikes. Thus, we study a model that is both realistic 
and general enough to enable us to incorporate different trends and extreme events. At the same time, 
as it will be seen in further sections, the model allows us to get closed form solutions for hedging and 
investment problems. 

Sections [2]l3] deal with the properties of such processes and the mathematical model of the market. Among 
the relevant results, we construct a unique martingale measure based on Girsanov's theorem. This measure 
guarantees the absence of arbitrage in our setting and shows that, under some scaling normalization, our 
model converges to that of Black-Scholes in distribution. In the next section, devoted to perfect hedging of 
options, we derive the fundamental equation, which, unlike the classical Black-Scholes model, is hyperbolic. 
In Section [5l we calculate the price of a call option and its simplifications including the case of Merton's 
model. In section [6l we study the questions relevant to imperfect hedging in the context of the considered 
model. 

Telegraph processes have been studied before in different probabilistic aspects (see, for instance, Gold- 



stein 



13i |, Kac [15|, Orsingher [19(1 and Zacks [2J]). These processes have been exploited for stochastic 



volatility modelling (Di Masi et al 3], Melnikov et al as well as for obtaining a "telegraph analog" of 
the Black-Scholes model (Di Crescenzo and Pellerey [6|). In contrast with the paper by Di Crescenzo and 
Pellerey, we use a more complicated and delicate construction of such a model to avoid arbitrage and to 
develop an adequate option pricing theory in this framework. 



2. Telegraph processes and auxiliary results 

Let (O, J^, P) be a complete probability space and a = o"(t), t > be a right-continuous Markov process, 
taking values ±1 and having intensities X± > 0: 

P{a{t + At) = +1 I a{t) = -1) = X-At + o(At), (1) 
P{a{t + At) = -1 I a{t) = +1) = X+At + o{At), (2) 

as At 0. The initial state a = (t(0) of the process (T{t),t > is deterministic and equal to +1 or —1. 

Further, we will consider all processes adapted to the filtration F = (F^)t>o (Fq = {0, 0}), generated by 
cr{t), t > 0, starting at fi, a = a{0) = ±1. We suppose that the filtration satisfies the "usual conditions". 
Recall (see e. g. Karatzas and Shreve [l^) that a filtration (F^) is said to satisfy the usual conditions 
if it is right-continuous and (Fq) contains all P-negligible sets of F. Under this assumption, the process 
c(t), t > can be viewed as a Markov flow controlled by transition probabilities (2.1)-(2.2). The time 
intervals tj — tj^i, j = 1, 2, ... (tq = 0), separated by instants of value changes tj = Tj, i = 1, 2, . . . 
are independent. Also, we denote Po- the conditional probability with respect to the initial state a, and E^- 
the expectation with respect to Po-. 

Let N"'{t) be the number of switches on [0, t] of the process cr{t), t > starting at o", a = a{0) = ±1. 
Note (see, for instance, Ross (23I ]) that N"' = N^{t), t > is a Poisson process with alternating intensities 
Xcn A_cr, A(j, . . ., 0" = ±1. 
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X^it) = J c^is)ds (3) 


and 

where cj(rj— ) is the left-limit of a{t) at Tj. As before, the superscript a = a{0) indicates the starting 
condition of the process cr{t), t > 0. 

Definition 2.1: The processes X'^ , a = ±1 are called (inhomogeneous) telegraph processes with states 
(c_, A_), (c+, A+) with starting at a. 

For A_ = A+ and — c_ = c+ = c, the process = ac J q (— ll^^'^^ds, t > is usually referred to as 
(integrated) telegraph process (see Goldstein [l^ and Kac |3])- 
The following lemma is evident. 

Lemma 2.2: Let the processes X'^ = X'^{t) and X^ = X"'{t), t > be telegraph processes with states 
[c±, X±) and {c±, X±) respectively, governed by the common Markov process a = cr{t). Then they are 
linearly connected X"'{t) = aX°'{t) + bt, t > 0, where 

C+-C- C+C- - C-C+ 

a = ac = , b = bc = , (5) 

C+ — C_ C4. — c_ 

and ac- + b = C-, ac+ + 6 = c+ . 

The next theorem could be considered as a version of the Doob-Meyer decomposition for telegraph 
processes with alternating intensities. 

Theorem 2.3: Let X^ be the telegraph process with states (c_, A_) and (c+, A+), defined in (0), and 
J'^ be the jump process, defined in Q), a = ±1. Then X^ + J" is a martingale if and only if 

A_/i_ = — c_, A+/i_|_ = — C4.. (6) 

Remark 1 : In particular, any (nontrivial) telegraph process without jumps (i.e. if h± = 0) never possess 
the martingale measure. 

Proof : In the particular case \± = A, h± = h, c± = c, the theorem follows from the martingale property 
of N(t) — At, t > 0. In a general situation, we have 

1 _ e-Mt-s) 

E{r{t) I F^) = r{s) + jH{t -s) + A,(,)a<,(,) , (7) 



1 _ e-^(*-*) 

E(X-(t) I F-)=X-(s)+5(t-s) + A,(,)d,(,) , a = ±l. 



(8) 



Here H = h^^h^^ A = A_ + , 7 = ^ , g = ^'^ , and ao 



Xa-ha-^^ — a-h.- 



-, da 



n — ' 



a = ±1. 
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To ensure the martingale property, we have to solve the equations 



7^^ + 5 = 0, 

a_ + d_ = 0, 
a+ + d+ = 0, 

which clearly leads us to dG]). □ 
Next, we study the properties of telegraph processes under a change of measure. Let be the telegraph 

N'it) 

process with the states (c^, A±), and = — Yl, (^*(j{t ~)I ^^{rj-) be the jump process with jump values 

i=i 

hj^ = — Cj_/A± > —1. Consider a probability measure P* with a local density with respect to Po-: 

Z'^{t) = '^\t=£t{Xl + j:), 0<t<T. (9) 

Here denotes the stochastic exponential (see e. g. Melnikov et al [l7|). 
Using properties of stochastic exponentials, we obtain 

Z-(t) = e^"W<(t), (10) 

where <(t) = H (1 + AJ,"(s)) with AJ-(s) = J^is) - J-(s-). 

s<t 

The process = Kj(t) can be represented as Kj(t) = Here a = ±1 indicates the initial 

direction, and the sequence k* ^ is defined as follows: 

<,<7 = <-l + " > 1> '^0,a = 1- (11) 

Thus if n = 2k, 

*\k / 



K*={i+h:ni + h*_^r, 



and if n = 2A; + 1, 



=(i+/i:)^+i(i+/.i. 



Theorem 2.4 : (Girsanov theorem) Under the probability measure P* , 

• the process = N^{t), t > is a Poisson process with intensities \*_ = \- — c*_ = A_(l + h*_) and 

a; = a+-c; = a+(i + /i;). 

• the process = X^[t), < t <T is a telegraph process with states (c_, Al) and (c+, A!j_). 

Proof: Let iTi^\t) = PaiN^'it) = n) and 7ri'^2(t) = P;(iV'^(t) = n), n = 0, 1, 2 where a indicates 
the initial state. The probabilities 7r^\t), a = ±1 solve the system 

1 vr^T^ \t=o= 0, n > 1; vr^'^^ \t=o= 1. 
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i:^:}{t) = E.(Z'^(t)l{^(,)=„}) = / e-+^V^r)(x, t)dx (12) 



-oo 



with a and b defined in ([5]), with c± replaced by c^. Here p^n \ n > are the (generalized) probability 
densities of the current position of the process X"{t)^ which has n turns with respect to the measure Po-- 
That is, for any borelian set A 

P^iX'^it) G A, N'^it) = n) = j P^:Kx, t)dx. (13) 

A 

Note that densities pn\x, t), cr = ±1 satisfy the equation 

^ + c.^ = -A.pM + A.pt^, n>l. (14) 
To prove the above, we note that, conditioning on a switch at the time interval (0, At), 

pI^\x, t + At) = (1 - X^At)pi"\x - c^At, t) + X^Atpl;_^-i{x, t) + o(At), At ^ 0. 
The system ()14p should be supplied with zero initial conditions p^^ \t=o= 0, n > 1 andpl^^ = e~^''*(5(x — 

Cat). 

Exploiting the equation (jl4p . we obtain from ()12p 



d 

^ = (ft - A. + acMiit) + A.(i - c:/xM;:liit)- 

The following evident equalities complete the proof: 

b - Xa + aca = cl - Xa = -A* , 
p(rMt=o= 0, n > 1, pi""^ \t=o=S{x). 

□ 

Remark 2 : System p4p could be solved exactly. Probability densities pi^^ , n > have a view pi^^ = 
g(-A±+!/t;±)t-ra^{^±)^ where u = AX/Ac = (A= — A__)/(c+ — c_) and are defined as follows (see ^^): 
q^^^ = 6{x — c±t) (and hence p^^^ = e~'^='=*5(x — c±t)), 

(+) _ K>^^- [c+t-xr-^jx -c.tr (-) _ A!^A!! (c+t-x)"(x-c„t)"-i 

(Ac)2"' (n-l)!n! *' '^^n (Ac)2n " n!(n-l)! *' ^ ' 



n = 1,2, 



and 



(±) A!^+^A!i {c+t-xr{x-c-tr i±) _Xl>y^ (c+t - x)"(x - c.t)" 

fcn+l - (Ac)2n+1 (^^1)2 - (Ac)2n+1 (^^1)2 ' 
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n = 0,1,2,... 

with Ot = Otix) = l{c_t<x<c+t}- 
This observation permits us to obtain the expUcit distribution p^^\x, t)dx = P±{xW(t) G dx}: 

oo 

= = e-^±* • 6{x - c±t) 

n=0 



+ e- 



-{\t+ux) 



y ^/ {c+t - x){x - C^t) 



(17) 



Here Io{z) = X^^q ^^/nl)^ ^^'^ ^li^) — ^oi^) modified Bessel functions, A = A± — i^c± — ^ '^"'"'^ 
This formula coincides with the main result of [2] (see Theorem 4.1). 



c+— c_ 



3. Market model based on telegraph processes 

Now we are ready to introduce the telegraph market model. The price of a risky asset S{t) follows 

dS{t) = Sit-)d {X''{t) + J^it)) , t>0 (18) 

and the process S{t), t > is right-continuous. Here we let X°' = X'^it), t > be a telegraph process 

with the states (c_, A_) and (c_|_, A+), c+ > c_, and J"" = J'^(t) = Yl ^o-(r -) with h± > —1. The 

i=i 

initial state of the market is defined by a = <t(0). 
The price of the non-risky asset (bank account) has the form 



i?(t)=e^'W, Y'^{t)= I r,(,)ds. 



(19) 



where r_, r+ > 0. 

According to the properties of stochastic exponentials, from (jlSp we have that 



S{t) = SoSt {X" + J") = Soe^'WK'^(t) 



(20) 



where 5*0 = 5(0) and 



s<t 



The sequence Kn,a, > is defined in (jlip (with h± instead of /ij_). 
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We assume that the parameters of the model (|18p -(|19 p satisfy the conditions 

^>0, a = ±l. (21) 

Under such conditions, we can find a unique martingale measure in the framework of the market (jlSp - 
(fT9]l . Recall that the measure P* , equivalent to Po-, is a martingale measure if the process {B{t)~^ S{t))t>Q 
is a P* -martingale. We define this measure by the density Z'^{t), t >0, ([9]), with hj_ = — Cj_/A±. 

Theorem 3.1 : Measure P* , defined by is the martingale measure if and only if 

C„ = Xa + 7 , CT = ±1. 

Moreover, under the probability measure P* , the process N'^ is a Poisson process with alternating intensities 

Proof : First of all, note that the condition (j21[) guarantees these inequalities: h% = — c* /Ao- = — 1 + {r^ — 
Ca)/{KK) > -1 and A* = Act - c* = (r^ - c„)/h„ > 0, a = ±1. Therefore, the process = Z''{t) = 
£t{X'^ + J") with c* = A(j + '^"h^" and /i* = — c*/Ao- defines the density of the new probability measure 
correctly. 

According to Theorem 12.41 the process X'^ — V is a telegraph process (with respect to P*) with the 
states {ca - ra, Xa - c*), a = ±1. From TheoremESl it follows that X'^(t) - ^ (t) + J'^it), t > is the 
P*-martingale if and only if 

(Act - Cl)ha- = -(Cct - Tct). 

Hence c* = Act + [c^ — ra-)/hf^ and h* = —c*/Xa = — 1 + {r^^ — Ca)/Xaha, and the Theorem is proved. □ 



To operate with assets B and S, we will exploit the notion of a trading strategy (portfolio) as a pair 
of two predictable processes ipt and ipt, t > 0. Here ipt and ipt are the amounts of the risky and the 
risk-free assets held in the portfolio at time t. The capital of the strategy is F^''^ = iptS{t) + il)tB{t). 
Such a strategy is self-financing if dF^'^ = iptdS{t) + iptdB{t), and admissible if F^"'"^ > 0, t > 0. Any 
self-financing strategy with a non-negative capital is called admissible. We will operate only with trading 
strategies which are admissible. 

A trading strategy is called an arbitrage strategy (at time T) if its initial capital is zero and Pct(F^''^ > 
0) > 0. It is well-known (see Delbaen and Schachermayer [S]) that the existence of a martingale measure 
guarantees that the market does not admit arbitrage. Hence, according to Theorem 13. 11 the market model 
([I8|)-([I9|) is arbitrage-free. 



Remark 1 : It is widely accepted that the telegraph process X'^ has a persistent character. Therefore, if 
J'^(t) = 0, the market has arbitrage opportunities (see the model considered in Di Crescenzo and Pellerey 
@]). The corresponding arbitrage strategy is described below. 

Assume r± = for simplicity. Take numbers A, B such that Sq < A < B < Soe'^+'^. Consider the 
following strategy: buy the risky asset at time ti = min{t £ [0, T] : S{t) = A}, and then sell it at 
time {2 = min{t E (ti, T] : S{t) = A or S{t) = B}. This strategy has no losses at time ti, because ti 
coincides with the switching time of X with zero probability. Hence the strategy creates a positive profit 
with positive probability P{S(t2) = B}. 



Now, let us discuss the convergence of (jlSp to the Black-Scholes model. First M.Kac [ij] noticed that 
the t eleg raph equation tends to the heat equation when c, A — > 00 and c^/A 1 (see also E.Orsingher 



19|], [20]). Moreover the homogeneous telegraph process X{-) converges in distribution, as c, A — > 00, and 
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(? /\ — > 1, to the standard Brownian motion w{-) in C([0, T]; (—00, 00)) (equipped with the sup-norm and 
the (T-algebra generated by the open subsets). See details in N. Ratanov [21[. 

The following theorem provides a similar connection between stock prices driven by geometric telegraph 
processes and geometric Brownian motion. Seeking for simplicity we consider the symmetric case A_ = 
A+, C- = a — c, c+ = a + c. 

Theorem 3.2 : Let A_ = A-|- = A ^ 00, c 00, 

cVA ^ vl a'/X ^ vl (22) 

Let ^ and 



a + XB/2^H, (23) 



where B = ln [(1 + /i-)(l + /i^ 
Then 



S{-)^Soe^p{vw{-)+fi-), (24) 



where means convergence in distribution in C([0, T]; (— cxd,oo)), and v = -s/v^+v^- 

Proof : Let f±{z, t) = E± '^'•'''Ht) \)q the moment-generating function of jump telegraph process 

xW(t) +lnKW(t), where n'^^^t) = Kn,± if A^W(t) = n (sequence i^n,±^ n > is defined by (llip with h± 
instead of h*^). We prove here the convergence 

/±(z, t) expifizt + v^zH/2), (25) 

which is sufficient for the convergence of one-dimensional distributions in ()24p . From Lemma [2. 2 1 it follows 

that 

00 „ 

f±{z, t) = E±e^K?'W+at+ln.c<±)W) ^ga.t^ / e-(-+ln.„,±)^(±)^^^ (26) 

where X^^"^ denotes a standard telegraph process with states (±1, A), and ph^\ n > here denote 
the (generalized) probability densities of x'^^\t): pi^^(x, t) = e~^^q^\x, t), where are defined in 
(dSD-dini) with c± = ±1 and A± = A. 
As /i± ^ 

f±(.z, t)r^e-^' / e--^e^-^/2pW(x, t)dx (27) 

Here and below f ^ g denotes the existence and equality of limits of / and g under scaling (|22 p - ()23p . 
Taking into account formulae (|15p - (jl6p and changing variables in the integral in (j26p . x' = cx, we obtain 



f±{z, t) ~ e'^"*+(^-^)* j e"^p(±)(x, t)dx, 
—00 

where p^^\x, t) is the density of the telegraph process X^^\t) with states (±c, A), A = Ae^-^/^. 



Next, note that 
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XzB Xz^B^ 

From (I22|)-(l23|), it follows that y/XB/2 ~ -a/VX and 

X — X + az ^ ^z + v'^z'^ /2. 



It is well-known that under the scaling c^/X the process X{t) converges in distribution to 

diffusion process v ■ w{t), and thus densities p^^\-, t) converge to the probability density of Vcw{t) (see 
e.g. i): 

Summarizing the above statements, we obtain ()25p . The convergence of any finite-dimensional distributions 



and the compactness property could be proved in the same manner as in N.Ratanov 2l[. □ 



Remark 2 : Condition ()23p in this theorem means that the total drift a + XB /2 is asymptotically finite. 
Here o = (c_ + Cj^)/2 is generated by the velocities of telegraph process X, and the summand XB/2 
represents the drift component (possibly with infinite asymptotics) that is provoked by jumps. If in (123p 
the limit of XB/2 is finite, then a ^ a = const, and in (124p the drift volatility term Va = 0. 

In general, by (f22]l - (f23]l . we have that ^fXB /2 — > —fa, so —VXB/2 represents the jump component of 
volatility. 



4. Fundamental equation and perfect hedging 

Consider a European option with maturity time T and payoff function f{S{T)). We assume / is a contin- 
uous and piecewise smooth function. To price these options, we need to study the function 



Fit, X, a) = E: e~'^^'-'^f{xe^'^^'-'^K%T-t)) , (28) 



(7 = ±1, < t < T, 

where E* denotes the expectation with respect to the martingale measure P*, under the initial state 
fj = fT(0) of the underlying Markov process o"(t), t > 0, a = ±1. 

Theorem 4.1 : Function F is a solution of the following hyperbolic system 

dF , , dF ^ 

— [t,x,a) + CaX—{t,x,a) 



= K + ^-^j-^)F{t, X, a) - ^-^-^F{t, x{l + h^), -a), 0<t<T, a = ±1 (29) 
with the terminal condition F{T, x, a) = f{x). 
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Proof: Note that V^it) = a^X'^(t) + with = K = "^l^J^'J^ (see LemmaE^]). Conditioning 

on the number of jumps, we can write 

oo „ 

F{t, X, a) = e-^--^^-*) ^ / e-''^yf{xeyKn,a)p^l{v, T - t)dy, a = ±1, (30) 

where n > are the (generahzed) probabihty densities of the telegraph process = X^{t), < t < 

T (with n turns), with respect to the martingale measure Pj_. The densities have a form ()13p - ()16p (with 
instead of A± and the same coefficient c±). Therefore the series in (|30p (and its derivatives) uniformly 
converges. It permits us to apply the equations (fl^ . Taking into account the identity Co-ar+br = r^, a = ±1 
(see Lemma l2.2p . we obtain from ()30p 

dF, , dF , 

-Q^{t, X, a) + CaX—{t, X, a) 

oo y 

= {r^ + X:)F{t, X, a) - A:e-^'-(^-*) ^ / e-'^^y f {xey KnM^-iiv ^ T-t)dy. 

By the equalities pT|) and A* = ^'f^'^' , the latter equation becomes ([29]) . □ 

Remark 1 : The system (I29p plays the same role for our model as the fundamental Black-Scholes equa- 
tion. In contrast with classical theory, this system is hyperbolic. In particular, it implies the finite velocity 
of propagation, which corresponds better to the intuitive understanding of financial markets. Note that 
the equations ([29|) do not depend on X±, just as the respective equation in the Black-Scholes model does 
not depend on the drift parameter. 

Now we consider the hedging problem for the option with a payoff function TC, which is FT-measurable. 
The self-financing strategy vr = {(ft, il^t), < t < T is called a hedge (perfect hedge, replicating strategy) 
if its terminal value is equal to the payoff of the option: 

F^ = n P-a.s. (31) 

For the wealth process Ft = F^ , we require that 

Ft = VtS{t) + iJtB{t), 0<t<T (32) 

and 

dFt = iptdS{t) + 4^tdB{t). (33) 
Let us rewrite ([33]) in the integral form 

t t jv-Ct) 

Ft = Fo+ / ^sS{s)dX%s)+ / i^sdB{s)+ ^ ^,^./i,(,^._)5(r,-). 



Using the equality -0* = B{t)^^{Ft — iptS{t)) (see the balance equation ([32]) ). we can rewrite the above 
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equation as 

t 

Ft = Fo + J r„^s)Fsds 



t N'it) 

+ J (PsS{s){c^(^s) -ra{s))ds + ^T,K{T,-)S{Tj-). (34) 

3=1 

To identify such a strategy in the case H = f{S{T)), note that Ft = F^ = B{t) E*[B{T)-'^n \ Ft] = 
F{t, S{t), cr{t)), where F{t, x, a) is defined by (j28|) and satisfies the fundamental equation ([29]) . Exploiting 
Ito's formula, generalized by Dolean and Meyer (see e. g. [13]), we get 

t t 
f dF r dF 

Ft = Fo+ I —{s,S{s),a{s))ds+ I —{s,S{s),a{s))S{s)c^^s)ds (35) 





Comparing the latter two equations and utilizing the fundamental equation (j29l) . we have (between 
jumps) 

Sit)c.it)^ + ^-r^is)F 

= 



F{t, S{t){l + K^t)), -a{t)) - F{t, S{t), a{t)) 



s{t)Kit) 

Moreover, from and ([35]) . we obtain the values of (prj'- 



(36) 



= P-.-P-.- = F{r„S{T,),a{T^)) - F{T„S{T,-),-a{T,)) ^ ^^^^ 



It turns out that the process (ft is left-continuous. To prove this, we note that based on (j20p . 

5(r,-)(l + /i.(.^._)) = 5(r,). (38) 
Now it is sufficient to apply §EI to (I36l)-(l371). 

5. Pricing call options 

The main goal of this section is to derive an exact formula for the initial price c of a call option with payoff 
(5(T) — K)~^ in the framework of the market ()18p - (jl9p . According to the theory on option pricing (see for 
example Karatzas and Shreve 16], Duffie @]), we have 



c'^ = E:{BiT)-'iSiT)-K)+), 



where K is the strike price and E* (•) is the expectation with respect to the martingale measure P*. 
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In case of the markovian model (|19p - (j20p . one can rewrite c*^ as 



c = = £ K{B{T)-\S{T) - K)+ l|;v^(T)=n}), (39) 



n=0 

where o" = ±1 indicates the initial state. 
We rewrite (1391) in the form 



c- = SoC/('^)(y,T)-K^x('^)(y,r) (40) 

with 



n=0 n=0 

where y = \\iK/Sq, h'^^ = InK^^o-, and functions uh \ Un \ n > are defined as follows: 

u'^^^y, t) = u'^^\y, t; A^, c±, r±) = E; [B{t)-H{x^^t)>y, N'{t)=n}] (41) 



e 



oo 

-b.t J e-«^-pg(x, t)dx 



y 



with Or = ^ and 6r- = — 



UtHy, t) = Ui-\y, t; Ai, c±, r±) = E; [B {t)-' StiX'^ + J'^)l{x^w>,, N'it)=n}] (42) 

CX) 

y 

The functions Un\y-, t), n> 1 satisfy the equations (see ([1] 



^{y, t) + c.^(y, t) = -(A: + r.)4'^)(y, t) + Kutliv, t) (43) 



with initial conditions Un \t=Q= 0, n > 1. As it is easy to see exploiting Remark [2] of Section 2 

4' 



these functions (un'^\ n > 1) are continuous and piece- wise continuously differentiable, and UQ^\y,t) 



|i=U— u, n ^ J 

.) are continuous 

e-(-^^+'"<T)t^^g^^ — y), o" = ±1. Moreover, \/n,u^^ = if y > c+t, and 

oo 

^^r^y,*) = pi'^HO = e-'"* / e-'^'-^pi':l(x, t)dx (44) 



if y < C-t. In the latter case, the system (|43|) takes the form 



^ = -(A: + r.)/>W+A:/>i-_f, n>l, (45) 
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pj,'") = e~{K+ra)t and p^n^ |t=o= 0, n > 1, ct = ±1. 

Lemma 5.1: The solution of system |^5[ j can he represented in the form 

p(r)(t) = e-(^*+-)*A(r)pM(t), CT = ±1, n>0, 
where A^T^ = (a;)[("+i)/2] (Al^)["/2] and functions are defined as follows: 



p(+) _ „-at 



p(-) = 1 



1+E 



k=l 



(n+l)fc ■ k\ 



, 0" = ±1, n > 1. 



(46) 



Here 



m, 



(+) = [n/2] , m(-) = [(n - l)/2] 



(m)fc = m{m + 1) ... (m + /c — 1), a = A!^ — Al + r+ — r_. 
Proof : Notice that in the particular case A!i_ = A!l = A and r± = the solution of system ()45p is well 



known: p^t\t) = 7Tn{t) = P{N{t) = n) = i^e 

Generally, we apply the following change of variables 

p(r)(t) = e-(^*+-)*A(r)pi'^)(t). 



(At)" „-At 



In these notations, we have P^'^\t) = e""*, a = {XX + r+) - {\*_ + r_)] ' {t) = I; P^ 
and the system 



(-), 



. p(±) 



|t=o= 



0, n > 1 



p{-) _ p{+) 



re > 1, 



(47) 



p(±) _ dP^±' 
~ dt ■ 

The latter system has the following solution 



2n+l — P2n+1 " (2n+l)! 



I _^ in+l)...{n+k) _ (-gf)'' 



p; 



{-) _ e 



2n — (2n)! 
>{+) _ i^" 



2n — ^2^. 



(2n+2)...(2n+fc+l) " fc! 

fc=l 

n(n+l)...(n+fc-l) (-gf) 
(2n+l)...{2n+fc) ' k\ 



1+E 

, I ^ (n+l)...(n+fc) (-gt) 

J- "T (2n+l)...(2n+fc) " fc! 
k=l 



which coincides with (H6 



□ 



Remark 1 : Formulas ()46p can be expressed using hypergeometric functions (Abramowitz and Stegun 



^^''^(i) = -1 ■ iFiim^n^ + 1; re + 1; -at), rez^ = [re/2] , rei^) = [(re - l)/2] 



re! 
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A hypergeometric function iFi{a; (3; z) is defined as 



^ , ^ . v-^ a(a + 1) . . . (q + n — 1) 
iFi(a; (3; z) = l + Y, .^^ 



1 + E 



(a)n 



^^n!/3(/3+ !)...(/? + n-1) ^ "!(/3)n 

Also using ([18]) . one can easily check that i-*2„^ ~ ^2n^ ~ o-P2n+i, n>0. 
Let us define the coefficients f3k,j, j < k: Pk,o = Pk,i = Pk,k-2 = Pk,k-i = 1, 



-z . 



Pk,j 



[i/2]! ' 



and the functions ipk,n- '■Po,n = -fbn+i and 



(48) 



k-l 



^k,n = ^a'-^~'(3k,,PL-j^ 1 < < 

j=0 



n. 



(49) 



For positive p, q, we define ^ = 0, Vq^' = e f^'^'' = Pi{p), cr = ±1, and for n > 1 



,(+) 



+1 = Hn+liP^ l) = P2n+l{p) + E |r'/'fc,n(P)' 



^2n 



(-) 



k=l 



^2n — ^2?i 
,(+) - 



{p, q) = P^Jip) + E fr'/'fc+i,n(p), 



n-l 



fe=l 



"2n - ^^2^^' Q) = P2n\p) + E frV^fc-l,™-! (p) • 

fc=l 



Now we can find expressions for ul^^ = u^\y, t) in the interval C-t < y < c+t. 
Theorem 5.2 : System admits a unique solution of the form 



(50) 



Wn\p, q), C-t <y < c+t, a = ±1, 
pl^\t),y <c-t, 



0, y > c+t, 



■u;/iere u;!'"^ = e-^^++'-+'>i~^^'~+'--'>PAl^\t\p, q), p = 
in Lemma\5.1\ 



(51) 



n > 0; functions pn'^ are expressed 



Proof: Evidently, ^^''^(y, t) = 0, if j? < 0, and u)a'{y, t) = pl^'it), if g < 0. For p, q > we have the 
system 



r dvi^^ _ J-) 



n > 1 



(52) 



with 
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and 



l.<o= 



e»^pW(p + g). (53) 



Here a = {\\ + r+) - {\*_ + r_) and Pi''^ n > 0, a = ±1 are defined in (gSD. 
It is straightforward to clieck tfiat tlie exact representation of the solution of (j52p for p, g > has the 

form ([50]) with (/?o,n = -P2n+i, V'l.n = -P^r^^ and 

^fc,n = V5fc-i,n-i, 1 < A: < n. (54) 

The proof is finished by the following proposition. 

Proposition 5.3: The solution of the system ( |5^| ) has the form 

fc-i 

j=0 



Proof Indeed, from (gg]) and (|47D it follows that 



k-l 
j=0 



By the identities P2n+i = -f'in+i ^'^'^ -^'in'* ~ -^in^ ~ a-P2n+i) > (see Remark [T] of this Section), we 
have 

j^O, j is even 



+ Yl a'-^-'PkAnlj-i- E a'-^Pk,jP2n-j. 

j>0, j is odd i>0, j is odd 

To complete the proof, it is sufficient to apply the identities Pk,2m+i = Pk-i,2m, Pk,2m — Pk,2m+i = 
Pk-i,2m-i, which are evident from the definition of Pk,n (see (H8l) ). □ 

Remark 2 : If A*_ = = A, r+ = r_ = r, then = tt^^ = 7r„ = ^e"^*, pi"^ = e-'-*7r„(t) and 
fk,n = Pta-k+i = (2n"-fc+i)! - Moreover, 



! — n \ / 



fc=0 

Remark 3 : It follows from (|5ip that functions Uq ^ and Uq^^ are discontinuous at g = and p = 

respectively. All other functions n > 1, defined in ()5ip . are continuous. The points of possible 

discontinuity of derivatives are concentrated on the lines p = and q = 0. For example, for and 
a = ±1, we have 



dq g=+o dq <i=-o 
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and 



y ^-{\%+r+)q _ 



dp p=+o dp p=-o 

Moreover, using (j5ip one can prove that G C"""*^. 

Similarly to system (H3]l for u"n \ the functions C/n*^^ = Un\y, t), n > 1, defined in (fl2|) . satisfy the 
equations 



+ c.^^ = -(A; + r,, - c.)[/W + A:(1 + K)U^-^^. (55) 



3t 9?/ 
Hence we obtain the following representation: 

\J^:\y, t; XI, c±, r±)=ul:Hy, t; A±, c±, 0), (56) 

where A^ = A* (1 + h„) = A* + - c^. 

Remark 4 : The formulas in (I40p have a different structure, which depends on the sign of ln(l + /i-)(l + 
h+). 

(i) If (1 + /i_)(l + h+) < 1, then ln(l + /i_) + ln(l + /i+) < and ft^T^ ^ -oo. The price of a cah option 
is given by the formula ()40p with 

and 

U = U^'^Hy, T) = u^^\y, T; A±, c±, 0), (57) 

where y = InK/So and n^'' = min |n : y — bj^^^ > c_t|, = min |n : y — bn^ > c+r| . 

(ii) If (1 + /i_)(l + h+) > 1, then ln(l + h_) + ln(l + /i+) > and ft^T^ ^ +oo. Denoting 



771^^ = max |n : y — h^"^ > c_t| , 

= max |n : y — 6^°^^ > c+t| , 
we obtain the call option price formula of the form ()40p with 



oo 



nW(y, T)= 5] 4'^)(y-6(-), T; Ai, c±, r±)+ P^^k^T), 

k=m[^^ A;=m'_f'+1 



and U^^^y, T) is defined in (l57D . 
Consider the following examples. 
Example 5.4 The Merton model. 
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= S{t-){cdt - hdN{t)), 

where N = N{t), t > is a (homogeneous) Poisson process with parameter A > 0. In this case, the formula 
(j40]) can be simphfied to 



c = SoU{lnK/So, T) - Ku{\nK/So, T). 

Here functions u and U are defined as follows. 

If < /i < 1 and c> r, then = bn = n ln(l — h) I — oo and 

no 

ti\nK/So, T) = e-^^ J^4^)(ln(K/5o) - 6„, T) 



u = u( 



(58) 



n=0 



e--^P(iV(T)<no)=e-'-^'vf„„(A*r), 



rT, 



where X* = {c — r)/h > and ^'^^(z) = e ^ In this case, the function U has the form 

n=0 



Uiy, T) = ^no{y{l-h)T). 



For h < and c < r, i. e. 51°^^ = nln(l — h) ^ +oo, we have 

n(y, T) = e-^^ (!-*„„ (A*r)) 



Uiy, T) = (A* (l-/i)r). 



In both cases. 



no 



inf{n : 5Qe"l'^(^-'^)+(^-'')^ > B{T)-^K} 



\n{K/So) - cT 
ln(l - h) 



Example 5.5 Let us consider another symmetric case A-|- = A_ := \, rj^ = := Cj^ = r + c,C- = r — c 
and /i+ = —h, h- = h; c > 0, < h < 1. These assumptions simplify the form of u^'^\ In this case we have 
A:^ = Al = c/h and 6^^^ ^ -oo. Here = nln(l - h^) and b^^^ = nln(l - /i^) +ln(l =F /i)- We denote 



ln(A75o) - {c±r)r 



ln(l - /i2) 

Function n*-^-* has the form (see Remark [2] of this Section) 



(59) 



u(+)(y,r)=e-(^/^+^)'^ 



2"- ^ TH/T x„ 2n+ [n/2] , 

^ n! ^ n! ^ U ' 

n=0 n=2n_+l fc=0 ^ 



where Pn = ^'^^^^'^^v^^"- > 0, qn = ^ ^" — ^ ^'^'^ ^ ~ lii(i^/S'o). Function ^ has a similar form. 
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6. Quantile hedging 



Consider an admissible strategy vr = (934, %l)t) with capital Ft = = (ptS{t)+iljtB{t) > for all t G [0, T]. 
For a given payoff function f(S{T)) and given admissible strategy {(ft, ipt) with the initial capital v, we 
define the set of successful hedging as 

A = A{v, ^, 4^, f ) = {ij : B{T)-'Ft > /}. 

If a strategy is a perfect hedge, then P„{A{c'^, ip, ip, /)) = 1, which requires the initial capital c'^ = 
E* [i?(T)~^/]. The problem of quantile hedging is to maximize probability of A under the budget restriction 



(60) 



P,^ {A{v, ip, tp, /)) max 

v<vo<E*, (i?(r)-V) = 

It is known (see Follmer and Leukert [lH) that ()60|) is equivalent to the following optimization problem 

Po-(^) max, 

< 

{B{T)-^f.lA)<vo. 



(61) 



Let A = Aa- be the solution of (j6ip . The perfect hedge {(p, ip) with initial capital vq for the claim 
/ = / • 1^ is the solution of ([60]) and its set of successful hedging A = A{v, (p, f) coincides with A. 
Moreover, the structure of the set A is 



A 



Using ©-([in]) and (P, we obtain 



dPo 

dp; 



T> 7 ■ / r ' 1 — const, 7 > 0. 



(62) 



dP 



dP. 



^ |t= £t{X^ + J:) = e^."(^)<(r) = e'^^'(^)+''^<(T), 



where a 



C-I-—C- 



and b 



c-^- — c_ 



Hence, the set of successful hedging A can be represented as 



For the standard call option with / = (5(T) — K)^ = {^Soe^' ^'^^ k!^ {T) — K^^ , the set A has the form 
A = i e-^"(^) > 7e^^<(r) [Soe^' ^%T) - [ = Q 



where 



An = |e-^"(^) > 7<,.e^^ (5o^n,.e^"(^) - k)^ , iV'^(r) = n} 
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An = {X%T)<yn, N%T)=n]. 

Here y„ = yn{l) = In Zn fen \ ^ — Ihk^^q-, and Zn — ^^(t) is the unique solution of the algebraic 
equation 

= 7<^,(K„,,)-"e^^(5oz - K)+. (63) 

It is clear that y„ = yn{l) decreases in 7 and y„ > IwK/Sq — hn'^ . 
To find the constant 7, we consider the equation 



Vo = SoY, [UIC\V - b^n \ T) - Ut\yn{l), T) 
n=0 
00 

-KY, [nW(y - h^:\ T) - nW(y„(7), T)] , (64) 



n=0 



where Un^ and C/n'^\ n > 0, cr = ±1 are defined in (|1T]) - (|12]) . y = IuK/Sq. As vq < c'^, the equation 
has a unique solution 7 = 7(^0) due to the monotonicity of y„ = yn{l)- Now the probability of maximal 
successful hedging set can be calculated as 

00 00 

P{A) = Y,^{An) = l-Y,^^nHyn{l), T; A±, c±, 0). (65) 

n=0 n=0 

Example 6.1 Let A4. = A_ = ^^^'^^ = ^^^^^ ■ It means that the initial measure P is the martingale 
measure, and the corresponding process is a homogeneous Poisson process with such intensity. Hence 
y„ = In ^'^^/^ — hn \ a = b = 0, and the equation ([Ml) for 7 = ^{vq) takes the form 



c{K, T) - ciK + 1/7, T) - -u^''\k + 1/7, T), 



■^0 

7 



where u^'^\z, T) = Yl Un — yn , T), and Un , n > are defined in ([5T]) . with Aj_ = X± and r-t = 0. 

n=0 

The probability of successful hedging is equal to 

p(i) = p^(5(r)<i^ + i/7) 



l-|;nW (in^^, T; A±, c±, o) , 



7 = 7(t;o). 



Let — a > 1, then we have 

A = {A-(r) < iV'^(r) = n} U {X'^(r) > N'^iT) = n] . 

Here yn^ = In z^^^ — b^^ and yl^"* = In Zn^ — bn'^ , where Zn^ and zi^^ are the solutions of ([63 
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The equation for 7 has the form 

00 

v = SoY^ [utHy - b^^\ T) - utHyi^\ T) + ^/('^^(yi^), T) 



n=0 



n=0 

and the solution of the quantile hedging problem is 

00 

P.(i) = l-^[4'^)(y«, T; A±, c±, 0)-4^)(yi2), T; A±, c±, 0)' 
The dual problem 



n=0 



v — > mm 



(66) 



(67) 



P<,(A(t;, ^, /)) > 1-e 

minimizes the initial capital under a fixed risk level. It can be solved as follows. Using ([65]) and (j66|) . we 
can find 7 from the equation Po-(j4'^) = 1 — e, i. e. 



Y.''n\yn{l). T; A±, c±, 0) =e (for - a < 1), 



(68) 



n=0 



J2[^nHy^n\l), T; A±, c±, 0)-n(r)(2/i2)(^), T; A±, c±, 0)' 

n=0 



(69) 



(for - a > 1), 

where y„ = In 2;„ — fe^T^ , and z„ = z„ (7) , n > solve the equation ([63]) . The set of successful hedging A is 
now defined and the optimal strategy is the perfect hedge of the claim / • 1^. 

Remark 1 : We would like to mention here a possible application of this type of hedging to risk- 
management of equity-linked life insurance contracts. The payoffs of such contracts depend on the evolution 
of risky assets as well as a survival status of insureds during the contract period [0,T]. Denote T{x) the 
remaining lifetime of a policy holder, who is currently of age x. Then the future payment may have the 
form max{ S'(T), ii'}/|T(3,)->T}., where K is the so-called maturity guarantee. It is natural to assume that 
T{x) does not depend on the stock market. Next, rewriting max{S'(T), K} = K + {S{T) — K)^, we reduce 
the pricing problem to that of a related call option. If = P{T{x) > T) is the survival probability, we 
can find the Brennan and Schwartz [5] price 



= E; [B{T)-^f . l{T(x)>T}] =. Px ■ E; [B{T)-^f] , 



(70) 



widely exploited in this area. This price is smaller than the fair price for the call option. Hence, perfect 
hedging is impossible with the help of underlying assets in the market. On the other hand, we can consider 
^Cx as vq in the context of quantile hedging above. For given T and x parameter ^p^ can be found from 
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actuarial life tables. Taking vq =^ Cx, we can construct the set A and the corresponding quantile hedge 
according to ([60]) - (f6T]) . 

Remark 2 : In the framework of the mod el (I18p -(|19 p it is possible to analyze optimal investment problems 
[l6l |. as well as shortfall risk minimization [l^. This analysis will be reported elsewhere later. Seeking for 
simplicity we omit these results here. 
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